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Abstract 

This  paper  presents  two  numerical  cases  of  disturbed  Navier-Stokes 
equations.  In  case  1,  four  order  accuracy  finite  difference  is  used  to 
investigate  the  evolution  of  negative  frequency  disturbance  in  Blasius 
boundary  layer,  results  show  us  that  negative  frequency  can  reduce  the 
streamwise  wavelength  of  T-S  wave,  and  accelerate  the  transition 
process.  The  effect  of  negative  frequency  lets  to  the  breakdown  of  fluid 
structure  becomes  more  significant.  And,  in  the  case  two,  two  order 
accuracy  finite  volume  method  is  used,  Reynolds  number  is  900,  and 
dimensionless  frequency  F is  86  in  this  case.  Results  present  that  the 
evolution  process  is  more  gradual  in  the  positive  frequency  compare  with 
the  negative  frequency  in  case  1.  And  the  pressure  contours  clearly 
present  a nonlinear  interaction  and  break  down  process  of  two  A 
vortexes. 


1.  Introduction, 

The  laminar-turbulent  transition  is  long  interested  by  researchers  of  fluid 
dynamics.  Stability  theory  has  been  developed  as  a theoretical  understand  of 
transition  phenomena,  its  early  period  is  the  stage  of  linear  theory,  L.  Prandtl1, 
W.  Tollmien  and  Grohne,  D 2,  W.  Heisenberg3  and  C.  C.  Lin4  greatly  contributed 
to  it.  The  character  of  this  period  is  the  development  of  theoretical  work  is  faster 
than  that  of  experiment  and  computation.  Up  to  1947,  by  the  using  of  ribbon 
technique,  G.  B.  Schubauer  and  H.  K.  Skramstad5  successfully  got  powerful 
experiment  supports  of  the  linear  stability  theory. 

Since  the  breakthrough  in  experiment  and  computation,  experiment  and 
computation  quickly  developed,  and  quickly  beyond  the  development  of 
theoretical  work,  because  the  difficulty  in  non-linear  analysis  obstructs  its 
advance.  The  most  important  achievements  in  the  stability  theory  during  this 
period  is  the  discovery  of  C-type6,  H-type7'8,  K-type9  and  N-type10, 11  breakdown. 
C-type  breakdown  can  be  explained  by  triad  resonance  theory,  and  H-type 
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breakdown  can  be  explained  by  the  secondary  instability  theory.  The 
experiments  of  Saric  and  Tomas12  show  that  for  the  different  amplitude  of  T-S 
wave,  fluid  will  transition  on  the  different  breakdown  route.  When  the  amplitude 
of  T-S  wave  range  from  0.2%  to  0.4%,  the  flows  will  present  C-type  instability, 
and  the  amplitude  of  T-S  wave  range  from  0.4%  to  0.6%,  H-type  instability  will 
appear  in  flow.  When  the  amplitude  of  T-S  wave  is  over  0.6%,  the  k-type 
instability  will  appear.  Since  the  secondary  instability  bases  on  the  assumption 
that  T-S  wave  is  independent  with  the  secondary  disturbance,  when  the 
amplitude  of  T-S  wave  is  high  enough,  this  condition  can  be  easily  satisfied.  This 
is  the  reason  H-type  breakdown  happen  on  the  higher  amplitude  T-S  wave 
conditions  compare  with  C-type  breakdown.  Since  the  clear,  directly  and 
universal  theoretical  analysis  about  Navier-Stokes  equation  has  been  unavailable 
up  to  now  because  of  the  nonlinear  difficulty,  numerical  simulation  become  the 
necessary  tool  to  the  stability  theory.  H.  Fasel13,  l4,  B.  J.  Bayly  & S.  A.  Orsag15 
and  C.  Liu  & Z.  Liu16  got  some  excellent  numerical  simulated  results  of 
transition  in  the  different  cases. 

The  purpose  of  this  paper  is  to  present  the  nonlinear  interaction  of  fluid 
structure,  and  effect  of  frequency  to  transition. 


2.  Governing  equation  and  numerical  method. 

Fluid  field  is  decomposed  into  two  parts:  mean  flow  and  disturbance.  The 
mean  flow  in  our  cases  is  Blasius  boundary  layer.  The  governing  equations 
become: 


M + j^duM+UtHufU^  + (dui(U2  + U°2)  + u?U2)  + + + 

3 t 3 £ dT]  dg 


...  3 3 3 . Am,  . 

*5e+,*5f+6a?J,-to-0 

3 u,  3m,  3m, 
dg  drj  dg 


(LI) 

(1.2) 


where  u,  are  disturbed  velocity  components,  p is  disturbed  pressure,  J is 
transformation  Jacobian,  uf  are  mean  flow  velocity  components.  Ui  and  Uf  are 

contravariant  disturbance  velocity  components  and  contravariant  mean  flow 
velocity  components,  and  hold: 

C/l  = — (Mi^  + M2^  + W3(^z ) , U2=~(u{fJx+U2Tly+UjTlz)  , 

U3=j(ulgx+u2gy+u3gz) 
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t/.°  = y («^V  + «2^v  + “3°£ ) . ^2  = y + «2\  + *#7*  ) 

^3  = y (MlVv  +1<2?v  +M3°?z) 

Inflow  displacement  thickness  is  used  as  dimensionless  constant.  The  following 
orthogonal  grids  are  used  for  the  numerical  simulation: 


y = 


1 

rimm°+y<mA'n™-'n) 


z = g (2) 


in  both  case,  yraax  = 40 , and  0 = 6.  And  the  grids  are  homogeneous  in  the 


streamwise  and  spanwise.  For  preventing  the  pressure-velocity  couple,  staggered 
grid  is  used. 

We  simulate  two  different  cases  for  the  description  the  evolutions  and 
interaction  of  T-S  wave  in  flat-plate  boundary  layer  transition.  The  first  case 
focus  on  the  effect  of  negative  frequency  disturbance  to  the  transition,  we  wish 
numerical  results  could  present  the  phenomena  that  the  negative  frequency 
disturbance  will  greatly  accelerate  the  breakdown  process.  And  for  this  purpose, 
finite  difference  method  with  second-order  backward  Euler  differences  in  time 
direction  and  fourth-order  central  differences  in  space  direction  is  used.  The 
inflow  is  the  superposition  of  solutions  of  O-S  equation,  and  holds: 

u.  =(j>(r])2idei(-m,)  +<f)(Ti)™ei(-(0'+p?)  +^(rp]d  + C.C.  (3.1) 

p = f (r7)2rf  el(-"n  + <t>\rj)3d  e^0,,+Pi)  + (j>\ri)3d  + C.C.  (3.2) 

The  side  condition  is  periodic,  which  hold: 

Qi.j,  1 = 0 i,j,NK+\  > 0 i.j.NK+2  = 01,/, 2 (4) 


So  every  periodic  in  spanwise  direction  will  be  investigated  by  K + 1 point.  And 
the  length  of  spanwise  is  set  as: 

_2n(NK  + l) 

n - 


PK 


(5) 


16 


For  avoiding  the  numerical  divergence  in  the  outflow  region,  the  previous  work 
use  buffer  zone  technique  in  the  outflow  region.  But  in  case  one,  we  replace  the 
buffer  zone  technique  with  the  first  order  upwind  scheme,  by  the  numerical 
dispassion  of  upwind  scheme,  the  numerical  divergence  is  successfully  avoided. 
The  using  of  upwind  scheme  bases  on  the  following  consider,  when  the  transition 
happen,  the  large  scale  vortex  structure  will  grow  up,  if  the  outflow  boundary 
locate  at  the  center  of  this  large  scale  vortex,  the  outer  flow  boundary  will  hold 
the  symmetric  boundary  condition,  such  as 


0/-1/2  0, 


/+1/2  ’ 


0 / —1/2  “ 0/* 


/+1/2 


where  / is  the  grid  of  outer  boundary.  In  this  condition,  we  have  the  following 
relation: 
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df  _ (/m -/m)-I^L./2 a^+i/2 +Hm,2 au,-1/2  214_i/2au,,i/2 

3£  A£  A£  A£ 

So 

ac/  3/  Bp  ay,  2K,„2ac/,.„,  ay, 

a/  a?  as  Rcaa?  a»  a?  as  Rea^ 

This  is  obvious  that  (6)  has  upwind  character;  this  kind  of  boundary  condition 
can  help  to  prevent  the  numerical  divergence. 

In  the  second  case,  the  disturbance; 

«,  =e12V(-a")  +efdea-m+Pi)  +efdei(-m-pi)  +C.C.  (7.1) 

u2=u3  = p = 0 (7.2) 

is  introduced  on  the  wall,  and  the  position  is  ( i',  1,  k ).  We  use  this  kind  of 
disturbance  to  simulate  the  behavior  of  disturbance,  which  is  excited  by  a ribbon 
on  the  wall.  The  two  order  spatial  accuracy  and  two  order  temporal  accuracy 
SIMPLE  method  is  used. 


2.  Results  and  Conclusions. 

In  the  case  1,  the  Reynolds  number  900  and  800  are  simulated,  and  P is  0.1. 
In  this  case,  we  set  (j)fd  =2 (j>fd  , the  wavelength  of  (3.1)  and  (3.2)  on  the 

. ■ P 

spanwise  direction  will  reduce  to  its  half,  — . Therefore,  a subharmonic,  which 

holds  a half  frequency  of  T-S  wave,  generates.  The  numerical  grid  of  case  1 is 
160x36x32.  Since  the  early  part  of  transition  is  a spatial  instability  dominant, 
wave  number  is  real,  so  we  assume  that  the  dimensionless  frequency  is  -86, 
according  to  analysis  of  A.  H.  Nayfeh  and  A.  Padhyen17,  the  group  speed  hold: 

— = r =^i| 
da  s dtW-s^e) 

where  X is  the  phase  angel  of  group  speed.  Therefore,  for  the  T-S  wave,  we  can 

get  the  following  formula: 

0(7?)ei(_fl>'+^)  =0(?7)exp[i(a^  -tujdr)]  = 0(77)exp[/(a<f;  (9) 

h to  L ge 


(8) 
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Our  numerical  results  in  case  1 (fig.l~fig.3)  proves  that  in  the  first  4 or  5 T-S 
wave,  the  group  speed  of  T-S  wave  is  homogenous.  Therefore,  from  (9),  it 
obtains  that: 

</>(ri)ei(~(0'+aS)  =</>07)exp  (10) 

where  stands  for  the  average  in  a streamwise  wavelength.  From  (10),  we  get 
a formula  of  streamwise  wavelength  lwave : 


/ 


U’tfl’tf 


= 2 7tl(are-(0 


(11) 


where  are  is  the  real  part  of  a . Although  the  wave  number  of  case  1 is  negative, 
the  propagation  direction  of  T-S  wave  and  the  direction  of  group  speed  are  still 
positive,  which  are  along  the  streamwise.  Therefore,  the  wavelength  lwave  of  case 
1 is  shorter  than  the  condition  of  corresponding  positive  frequency,  such  as  case 
2. 

Figure  1 , 2 and  3 are  consequences  of  case  1 , Figure  2 clearly  present  that,  in 
the  second  periodic  of  T-S  wave,  it  breaks  into  two  parts  in  the  spanwise,  this  is  a 
typical  harmonica  behavior.  As  the  evolution  of  T-S  wave,  structures 
continuously  break  down;  fluid  contour  becomes  more  and  more  irregular  and 
complex.  Turbulence  happens.  The  strength  of  the  amplitude  of  disturbance  in 

figure  3 is  higher  than  that  of  figure  1 and  figure  2,  and  max(02d ) = 0.04  in 

figure  3,  and  max(02d ) = 0.02  in  figure  1 and  figure  2,  where  <j)fd  is  the  profile 
of  coming  flow. 

Figure  4 to  1 1 are  numerical  result  of  case  2.  In  this  case,  Reynolds  number  is 
900,  f5  is  0.1,  and  the  dimensionless  frequency  F is  86.  But,  in  case  2, 

(j)fd  = 3(j>:d . The  numerical  grid  is  140x32x32.  Numerical  ribbon  sets  on  the 
grids  ( 30, 1,  k ).  The  maximum  receptive  velocity  on  the  ( 30,  2,  k ) is  0.005. 
Results  clearly  show  us  an  entire  interaction  process  of  two  X vortexes.  At  first, 
each  X vortex  will  be  stretched  by  the  other,  their  head  will  move  closely  in 
their  propagation  direction.  In  the  region  Re  = 1 200 , two  X vortexes  merge 
together.  And,  quickly,  breakdown  happens. 
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Figure  1,  Spanwise  velocity  contour,  Re  = 900,  F = 86.  y = 2.0 
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Figure  2,  Spanwise  velocity  contour,  Re  = 900,  F = 86.  y = 1 .4 


S tream  w ise 

Figure  3,  Spanwise  velocity  contour,  Re  = 800,  F = 86.  y = 2.0 
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Figure  4,  Pressure  contour.  Re  = 900,  F = 86.  y = 2.0,  t = 7T, 
Ribbon  position  x = 303.94 


Figure  5,  Pressure  contour,  Re  = 900,  F = 86.  y = 2.0,  t = 7.125T, 
Ribbon  position  x = 303.94 


Figure  6,  Pressure  contour,  Re  = 900,  F = 86.  y = 2.0,  t = 7.  25T, 
Ribbon  position  x = 303.94 
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Figure  7,  Pressure  contour,  Re  = 900,  F = 86.  y = 2.0,  t = 7.  375T, 
Ribbon  position  x = 303.94 


Figure  8,  Pressure  contour,  Re  = 900,  F = 86.  y = 2.0,  t = 7.  5T, 
Ribbon  position  x = 303.94 


Figure  9,  Pressure  contour,  Re  = 900,  F = 86.  y = 2.0,  t = 7.  625T, 
Ribbon  position  x = 303.94 


